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ABSTRACT 

This paper presents contour maps of the global geoid 
and free air gravity anomalies of the GEM-4 Earth 
gravitational model derived at the Goddard Space 
Flight Center in 1972. 
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THE GEOID AND FREE AIR GRAVITY ANOMALIES CORRESPONDING 
TO THE GEM-4 EARTH GRAVITATIONAL MODEL 


1.0 INTRODUCTION 

This paper presents the geoid and free air gravity anomalies corresponding to 
the GSFC GEM-4 model of the Earth's gravitational field. The GEM-4 solution 
(Lerch, et al. 1972) is based upon a combination of electronic, laser and optical 
observations of 27 satellites with 5° x 5° mean free air gravity anomalies cover- 
ing about 70% of the world. 


2.0 COMPUTATION OF THE GEOID 

The procedure for computation of the geoid consisted of first fixing a value of the 
potential of the geoid, W 0 , and then calculating the geoid height as 


N = r - r ( : 


( 1 ) 


where: 

r is the radial distance to the equipotential surface (the geoid) defined 
by W 0 and the potential coefficients of the GEM-4 gravity model. 

r E is the radial distance to a selected reference ellipsoid defined by a 
semimajor axis (a e ) and flattening (f). 

The radial distance, r, to the equipotential surface W 0 at a particular latitude 
and longitude 0, \ is determined by using equation 2 . 


W () = 0 (r , 0, X) = 


GM 

r 

oo 

"L 

£ (€ < 

C nm cos m\ + S nm sin mX) P nm (sin0) 


n=2 

m=o 

— 


to 2 r 2 


COS 2 0 


( 2 ) 


1 



where 


GM 

a e 


= the product of the gravitational constant and the mass of the 
earth 

= semimajor axis of the reference ellipsoid 


r 


= geocentric radius 


to 


C nm and S nm 
■^nm (sin#) 


earth's angular velocity 

fully normalized spherical harmonic coefficients of the 
gravitational potential 

the Associated Legendre Polynomial (fully normalized) 


The only unknown in this equation is r. The value of r at each specified point 
of latitude and longitude is obtained by iteration. Geoid undulations were then 
computed as the separation between the geoid and the reference ellipsoid by 
means of equation (1). 


The parameter values adopted for these calculations were: 

GM = 3.986009 x 10 s km 3 /sec 2 

a e = 6378.142 km 

to = 0.72921151467 x 10^* rad/sec 

l/f = 298.255 

W 0 = 6263687. 5 kgal m 


A map of the GEM-4 geoid heights plotted with a 2 meter contour interval is 
presented at the end of this report. 

3.0 COMPUTATION OF FREE AIR GRAVITY ANOMALIES 

A plot of the free air gravity anomalies (Ag) is also presented at the end of the 
report. These were obtained by evaluating the following equation. 


Ag = 7 (n-1) [c nm cos ml + S nm sin ml] P nm (sin#) (3) 


n=2 m=0 


2 



where 


7 = Mean value of gravity over the earth (0.98 x 10 6 mgal) 


In equation (3), the C 2 o and C 40 terms do not represent the complete coefficients 
but rather the difference between the complete coefficients and the coefficients 
compatible with the ellipsoid used in computing N . The difference values used 
wei'e AC 20 = 0.01954 x 10* 6 andAC 40 = -0.24 1 7 x 10' 6 (fully normalized). 
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